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Abstract. This paper ia a written version of a one hour lecture given 
on Lov Grover's quantum database search algorithm. It is based on , 
§, and §. 
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1. Problem definition 



We consider the problem of searching an unstructured database of N = 2 n 
records for exactly one record which has been specifically marked. This can 
be rephrased in mathematical terms as an oracle problem as follows: 
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2 SAMUEL J. LOMONACO, JR. 

Label the records of the database with the integers 

0,1,2, ... ,N-1 , 

and denote the label of the unknown marked record by xq. We are given 
an oracle which computes the n bit binary function 

/: {0,1}" — {0,1} 

defined by 

{1 if x = xo 
otherwise 

We remind the readers that, as a standard oracle idealization, we have no 
access to the internal workings of the function /. It operates simply as a 
blackbox function, which we can query as many times as we like. But with 
each such a query comes an associated computational cost. 

Search Problem for an Unstructured Database. Find the record 
labeled as xq with the minimum amount of computational work, i.e., with 
the minimum number of queries of the oracle f. 

From probability theory, we know that if we examine k records, i.e., if we 
compute the oracle / for k randomly chosen records, then the probability of 
finding the record labeled as xq is k/N. Hence, on a classical computer it 
takes O(N) = 0(2 n ) queries to find the record labeled xq. 



2. The quantum mechanical perspective 



However, as Lov Grover so astutely observed, on a quantum computer the 
search of an unstructured database can be accomplished in 0(y/N) steps, or 
more precisely, with the application of 0(y/NlgN) sufficiently local unitary 
transformations. Although this is not exponentially faster, it is a significant 
speedup. 



Let H2 be a 2 dimensional Hilbert space with orthonormal basis 

{|0>,|i» ; 

and let 

{|0>,|1>, ... ,\N-1)} 
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denote the induced orthonormal basis of the Hilbert space 

N-l 

H = 0H a • 



From the quantum mechanical perspective, the oracle function / is given 
as a blackbox unitary transformation Uf, i.e., by 

H®H 2 H®H 2 

\x)®\y) i — * \x) <g> \f(x) © y) 
where '©' denotes exclusive 'OR', i.e., addition modulo 2.[] 



Instead of Uf, we will use the computationally equivalent unitary trans- 
formation 



I {xo) (|x» = (-!)'<*> \x) 



\xq) if X = Xq 



\x) otherwise 



That I\ XQ } is computationally equivalent to Uf follows from the easily veri- 
fiable fact that 

and also from the fact that Uf can be constructed from a controlled- 



and two one qubit Hadamard transforms. (For details, please refer to [flOfl , 

go 



The unitary transformation L XQ \ is actually an inversion Q in TL about 
the hyperplane perpendicular to \xq). This becomes evident when L XQ \ is 
rewritten in the form 

I\x. ) = I ~2\x ) (x \ , 

where 'J' denotes the identity transformation. More generally, for any unit 
length ket \tp), the unitary transformation 

I m =I-2\4>) (V| 

is an inversion in TL about the hyperplane orthogonal to 

1 Please note that Uf — (y o t) (/), as denned in sections 10.3 and 10.4 of ]H| . 
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3. Properties of the inversion 



We digress for a moment to discuss the properties of the unitary trans- 
formation Ji^\. To do so, we need the following definition. 

Definition 1. Let and \x) be two kets in H for which the bracket product 
(ip | x) i s a rea l number. We define 

S c = Span c (|V> , |x» = {« IV>) + P |x) E H | a, (3 E C} 

as the sub-Hilbert space of H spanned by and |x)- We associate with 
the Hilbert space Sc a real inner product space lying in Sc defined by 

S R = Span R , |x» = {a + b \ X ) E H \ a, b G R} , 

where the inner product on 5r is that induced by the bracket product on 7i . 
If\ip) and \x) are also linearly independent, then 5r is a 2 dimensional real 
inner product space (i.e., the 2 dimensional Euclidean plane) lying inside of 
the complex 2 dimensional space Sc- 



Proposition 1. Let and \x) be two linearly independent unit length 
kets in 7i with real bracket product; and let Sc = Spanc , |x)) an d 5r = 
Span R (\ip) , |x)). Then 

1) Both Sc and 5r are invariant under the transformations I\ x \, and 
hence Lm o L x \ } i.e., 



I\i/>) ( 5 c) = 


Sc 


and 


I\i/>) (<5m) = <5k 


J lx> ( 5 c) = 


Sc 


and 


I\ x ) = s ^ 


I w I\ x) (S c 


)=S C 


and 





2) If L^x^ is the line in the plane 5k which passes through the origin and 

which is perpendicular to \ then 1^ restricted to 5r is a reflection in 
(i.e., a Mobius inversion about) the line Luiv. A similar statement 

can be made in regard to \x)- 

3) // l^" 1 ) is a unit length vector in 5k perpendicular to \ip), then 

- T M = T \^) ■ 
(Hence, (ip 1 - | x) is real.) 
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Proposition 2. // is a unit length ket in 7i, and if U is a unitary 
transformation on 7i, then 



4. The method in Lov's "madness" 



Let H : TC — ► TC be the Hadamard transform, i.e., 

n-l 

H = (g)tf( 2 ) , 
o 

where 




with respect to the basis |0), |1). 

We begin by using the Hadamard transform H to construct a state 
\ipo) which is an equal superposition of all the standard basis states |0), 
|1),. . . ,\N — 1) (including the unknown state \xq)), i.e., 

JV-l 

|lfo> = # |0> = -= J» • 

Both jV'o) and the unknown state \xq) lie in the Euclidean plane 5r = 
Spanu (\ipo) , \xq)). Our strategy is to rotate within the plane Sr the state 
|^o) about the origin until it is as close as possible to \xq). Then a measure- 
ment with respect to the standard basis of the state resulting from rotating 
|^o)) wm produce \xq) with high probability. 

To achieve this objective, we use the oracle I\ Xo ) to construct the unitary 
transformation 

Q = -HI\q)H~ 1 I\ Xo ) , 

which by proposition 2 above, can be reexpressed as 

Q = ~ / IV'o> I No> • 

Let \xq) and \4>q) denote unit length vectors in <Sr perpendicular to \x ) 
and jV'o), respectively. There are two possible choices for each of \xq) 
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and h/'o / respectively. To remove this minor, but nonetheless annoying, 
ambiguity, we select \xq) and \ipQ } so that the orientation of the plane 5r 
induced by the ordered spanning vectors IV'o); \xq) is the same orientation 
as that induced by each of the ordered bases \xq ), \xq) and \ipo), "0o )■ 
(Please refer to Figure 2.) 

Remark 1. The removal of the above ambiguities is really not essential. 
However, it does simplify the exposition given below. 









•l x o> 


a jt \ 




V„>/ 


— " > 







L IV 





Figure 2. The linear transformation Q| c is reflection in the line L\ ±\ 
followed by reflection in the line £u ) which is the same as rotation by the 
angle 2(3. Thus, Q\s„ rotates | V^o) by the angle 2(3 toward \xq). 

We proceed by noting that, by the above proposition 1, the plane 5r 
lying in 7i is invariant under the linear transformation Q, and that, when 
Q is restricted to the plane 5k, it can be written as the composition of two 
inversions, i.e., 

Q\s R = I \i,^) I \x ) ■ 

In particular, Q\ s is the composition of two inversions in 5r, the first in 

the line L\ ±\ in 5r passing through the origin having \xq) as normal, the 

ro / 

second in the line L\^p ) through the origin having IV'o") as normal.^ 
We can now apply the following theorem from plane geometry: 



2 The line ^U-M is the intersection of the plane Sm with the hyperplane in Tt orthogonal 
to \xo). A similar statement can be made in regard to L\^ a ). 
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Theorem 1. If L\ and Li are lines in the Euclidean plane M 2 intersecting 
at a point O; and if (3 is the angle in the plane from L\ to L<i, then the 
operation of reflection in L\ followed by reflection in Li is just rotation by 
angle 2(3 about the point O. 



Let (3 denote the angle in Sk from Li±\ to £|^ ), which by plane geometry 
is the same as the angle from |xg~) to \ipo), which in turn is the same as the 
angle from \xo) to IV'o")- Then by the above theorem Q\ Sr = I^±^I\ Xo ) is 
a rotation about the origin by the angle 2(3. 



The key idea in Grover's algorithm is to move IV'o) toward the unknown 
state \xq) by successively applying the rotation Q to IV'o) to rotate it around 
to \xo). This process is called amplitude amplification. Once this 
process is completed, the measurement of the resulting state (with respect 
to the standard basis) will, with high probability, yield the unknown state 
|xq). This is the essence of Grover's algorithm. 



But how many times K should we apply the rotation Q to IV'o)? If we 
applied Q too many or too few times, we would over- or undershoot our 
target state \xq). 



We determine the integer K as follows: 



Since 

IV'o) = sin (3 \x ) + cos (3 Xq / 
the state resulting after k applications of Q is 

iV'fc) = Q k IV'o) = sin [(2k + 1) (3] \x ) + cos [(2k + 1) (3] 
Thus, we seek to find the smallest positive integer K = k such that 

sin [(2k + 1) (3] 

is as close as possible to 1. In other words, we seek to find the smallest 
positive integer K = k such that 



(2k + l)(3 
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is as close as possible to tt/2. It follows that| 



vr 1 

K = k = round I — — 



4f3 2 / 

where "round" is the function that rounds to the nearest integer. 



We can determine the angle (3 by noting that the angle a from \ipo) an d 
|xo) is complementary to 0, i.e., 

a + f3 = tt/2 , 

and hence, 

1 7T 

(xq I ifio) = cos a = cos(— — p) = sin/3 . 



Thus, the angle (5 is given by 



13 = sin 1 ( | — L= (for large iV) 



and hence, 



K = k = round ( -, r I ~ round ( — \/iV ] (for large N). 



5. Summary of Grover's algorithm 



In summary, we provide the following outline of Grover's algorithm: 
3 The reader may prefer to use the floor function instead of the round function. 



GROVER'S ALGORITHM 



9 



Grover's Algorithm 



STEP 0. 



STEP 1. 



STEP 2. 



(Initialization) 

N-l 
3=0 

k < — 



Loop until k = round 



4 sin 



k 



-HI\o)HI\ Xo) 



round (jVN -\ 



Measure \ip) with respect to the standard basis 

|0) , 1 1) , ... , | N — 1) to obtain the marked unknown 

state \xq) with probability > 1 — ^. 



We complete our summary with the following theorem: 



Theorem 2. With a probability of erroif\ 

Prob E < —, 

Grover's algorithm finds the unknown state \xq) at a computational cost of 

O ( VnIrN^ 

Proof. 



Part 1. The probability of error ProbE of finding the hidden state \xo) is given 
by 

Prob E = cos 2 {(2K + 1)13} , 



where 



K 



sin 



l / i 



round f ^ — \ 



If the reader prefers to use the floor function rather than the round function, then 

A. — 4 



probability of error becomes ProbE < 4 
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where "round" is the function that rounds to the nearest integer. 
Hence, 

^-1<K<^ => Z-p<(2K + l)0<* + (3 

=> sin (3 = cos (§-/?)> cos [{2K + 1) 0\ > cos (f + (3) = - sin/3 

Thus, 

Prob E = cos 2 [(2if + !)/?]< sin 2 (3 = sin 2 ^sin' 1 (^=j ) = ^ 



Part 2. The computational cost of the Hadamard transform H = 0q _1 -^ 2 ' ) 
is O(n) = 0(lgN) single qubit operations. The transformations —I\o) 
and I\xq) each carry a computational cost of 0(1). 

STEP 1 is the computationally dominant step. In STEP 1 there are 
O (^/N^j iterations. In each iteration, the Hadamard transform is ap- 
plied twice. The transformations — and I\ XQ \ are each applied once. 
Hence, each iteration comes with a computational cost of O (lgTV), and 
so the total cost of STEP 1 is O(^NlgN). 

□ 



6. An example of Grover's algorithm 



As an example, we search a database consisting of N = 2 n = 8 records 
for an unknown record with the unknown label xq = 5. The calculations for 
this example were made with OpenQuacks, which is an open source quantum 
simulator Maple package developed at UMBC and publically available. 



We are given a blackbox computing device 



In 



'|?> 



Out 
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that implements as an oracle the unknown unitary transformation 



l \*o) 



'|5> 



/ 1 

10 

10 

1 







y o o o o 



\ 







-10 

10 

10 

1 / 



We cannot open up the blackbox 



f|?> 



to see what is inside. 



So we do not know what I\ XQ ^ and xq are. The only way that we can glean 
some information about xq is to apply some chosen state \tp) as input, and 
then make use of the resulting output. 



Using of the blackbox 



a computing device 
operator 



- HI \o) HI \?) 



as a component device, we construct 
-> which implements the unitary 



Q - - HI \o)HI\ xo ) - 



( 


-3 


1 


1 


1 


-1 


1 


1 


1 \ 




1 


-3 


1 


1 


-1 


1 


1 


1 




1 


1 


-3 


1 


-1 


1 


1 


1 




1 


1 


1 


-3 


-1 


1 


1 


1 




1 


1 


1 


1 


3 


1 


1 


1 




1 


1 


1 


1 


-1 


-3 


1 


1 




1 


1 


1 


1 


-1 


1 


-3 


1 


V 


1 


1 


1 


1 


-1 


1 


1 


"3/ 



We do not know what un itary transformation Q is i mplement ed by the 
device - 



-HI\q)HI\ ? ) 



because the blackbox 



<|?> 



is one 



of its essential components. 



STEP 0. We begin by preparing the known state 



|V>o> = H\0) = (1,1, 1,1, 1,1, 1,1) 



transpose 
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STEP 1. 1 We proceed to loop 

K = round 



TT 



4 sin- 1 (1/V8) 2 



times in STEP 1. 
Iteration 1. On the first iteration, we obtain the unknown state 



|Vi> = Wo) = (M, 1,1,5,1,1,1) 



transpose 



Iteration 2. On the second iteration, we obtain the unknown state 



l^2> = Q IV>i) = ^ (-1, -i, -i, -i, 11, -i, -i, -i)*™*°* e 



and branch to STEP 2. 



STEP 2. We measure the unknown state \ip2) to obtain either 

|5) 

with probability 

121 

Probsuccess = sin 2 {{2K + 1) /?) = —- = 0.9453 
or some other state with probability 

Prob Fa uure = cos 2 {{2K + 1) p) = ^ = 0.0547 
and then exit. 
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